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We investigate theoretically the phenomenon of so-called fast light in an unconventional regime, using pulses 
sufficiently short that relaxation effects in a gain medium can be ignored completely. We show that previously 
recognized gain instabilities, including superfluorescence, can be tolerated in achieving a pulse peak advance of 
one full peak width. © 2008 Optical Society of America 
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1. Introduction 

Studies of electromagnetically induced trans- 
parency (EIT) 1 have suggested ways that a strong 
laser field can dress a dielectric medium which can 
result in a modification of the absorptive and dis- 
persive characteristics of the medium. This alters 
the group index of the medium such that group ve- 
locity of a weak probe pulse injected in the dressed 
medium can be modified. 2 4 Observations of slow 
light 5 and fast light 6 in such dressed media have 
been reported. This modification of the medium, 
and subsequent change in the group velocity leads 
to advanced or retarded propagation of a pulse 
peak. Reviews of early fast light experiments can 
be found in Refs. 2,3,7. 

In this paper we are concerned with fast light, 
which can occur if the group velocity of a pulse 
is greater than the vacuum speed of light, or even 
negative. High dispersion allows for large change 
in the group index, which leads to strongly altered 
group velocity. One particularly effective technique 
in obtaining high dispersion is to manipulate the 
spectral domain between two resonance lines with 
external cw fields. This can create a narrow spec- 
tral window with high dispersion, through which a 
long probe pulse will propagate with strongly al- 
tered group velocity and little or no loss. A steady- 
state high dispersion window can be created in 
more than one way, 8 but the narrow bandwidth 
of the achieveable spectral window forces narrow- 
bandwidth restrictions on the probe pulse. Noise 
limitations have also been discussed. 9, 10 

Coherent propagation of very short pulses is not 
affected by these narrow-line considerations. We 
will treat pulses short enough that upper level re- 
laxation can be ignored and steady state cannot be 



reached. Previous theoretical examinations of such 
fast-light short pulses in gain media have been re- 
ported. 11 ' 12 However because of instabilities asso- 
ciated with gain media, large pulse advance was 
not thought possible. 

We have recently theoretically demonstrated 13 
that sech-shaped pulses can give a substantial 
pulse advance of more than 10 pulse widths even 
when including instabilities associated with the 
two- level Area theorem. 14 This treatment did not 
include the quantum instability of superfluores- 
cence, first predicted by Dicke in 1954, 15 and ob- 
served in the late 1970's. 16 ' 17 Superfluorescence 
theory predicts that a group of atoms will spon- 
taneously decay collectively at a rate much greater 
than the single atom's decay rate. This enhanced 
spontaneous emission rate causes the coherent gain 
medium to relax very quickly to the ground state, 
potentially greatly interfering with superluminal 
pulse propagation. Here we describe the results of 
a detailed examination of fast light in a fully co- 
herent gain medium and show that a sech-shaped 
input pulse still confers advantages leading to a 
significant pulse advance with both types of insta- 
bility accounted for. 

2. Theoretical Model 

To keep our model as simple as possible, while 
retaining the physics of a fully coherent inverted 
medium, we assume the medium consists of a col- 
lection of inverted two-level atoms confined to a 
cell of length L and cross sectional area a as shown 
in Fig. [H Our approach is similar to that taken 
by earlier workers. 11, 12 We allow a range of detun- 
ings A, so the atoms are subject to inhomogeneous 
broadening. The Hamiltonian of a single atom in 
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Fig. 1. Two- level atom with level 1 connected to 
level 2 via the Rabi frequency Q of a laser field de- 
tuned from resonance by an amount A, and sketch 
of significant pulse advance in an inverted two-level 
gain medium of length L, and cross sectional area 
a. 



inhomogeneously broadened Beers absorption co- 
efficient, and g = N{i 2 uj/coh. Eqns. ([2]) and ([3]) 
combine to form the well known nonlinearly cou- 
pled Maxwell- Schrodinger (MS) system. 

One can see the differences between EIT theory 
and this description. We are assuming a completely 
coherent interaction between the laser field and the 
atoms by neglecting all relaxation terms in Eq. ([2]) . 
This approach is valid when the temporal width of 
the probe pulse and the pulse transit time is short 
enough to ignore relaxation. 



the rotating wave picture is given by 



F=-^|l)(2|-^|2)(l| 
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where A = uj 2 i — w is the detuning of the laser 
below resonance, UO21 is the frequency difference 
of levels 1 and 2, Q = 2/jlS/H is the Rabi fre- 
quency associated with the 1 — ► 2 transition, \i 
is the dipole matrix element, £ is the complex en- 
velope function of the electric field. The electric 
field has been decomposed into a slowly varying 
envelope function and a carrier frequency such that 
E = £ e - l ( kx -ut) _|_ c c>5 where uj is the laser fre- 
quency, and k = uo/c is the wavenumber. For sim- 
plicity we will restrict our attention to Rabi fre- 
quencies that are purely real for the remainder of 
the paper. 

We assume a fully coherent atomic state of the 
form \i/j(x, t)) = ci(x, i)\l)+C2(x, t)\2) with c\ and 
C2 being the corresponding probability amplitudes. 
Schrodinger's equation along with the Hamiltonian 
in Eq. (pQ) then gives 
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Maxwell's wave equation written in the slowly 
varying envelope approximation (SVEA) yields an 
equation for the pulse envelope via the Rabi fre- 
quency: 



c— + — V 

dx dt , 



where F(A) = ^e - ^ A )"/ 2 , and T 2 * is the inho- 
mogeneous lifetime. The product 7rgF(0)/c is the 
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3. Analytic sech Pulse Solutions 

We use a mathematical approach to the coupled 
MS system that provides analytic expressions for 
fast-light pulses. The method we use, a variant of 
the Backlund transformation technique developed 
by Park and Shin, 18 is flexible enough that we can 
use it to treat different possible experimental situ- 
ations such as segmented and stacked media, 19 as 
well as media that have a variable density of atoms. 

To model a medium confined to an atomic va- 
por cell between positions xq and x\ with length 
L = x\ — xq as shown in Fig. [TJ we take the 
simplest special case for a variable density, i.e., 
atoms with density TV inside the medium, and 
N = outside the medium. We will assume that 
the atoms in the medium are uniformly excited: 
\ip(x,t = —00)) = |2), a state of perfect inver- 
sion. The exact SVEA pulse solutions are given by 
the following segmented expression that describes 
a fully continuous pulse in time and space: 



{fsech^-^) x<x 
^sech^ (t - f- g + O ) x < x < x x 
^sech^ (t- z + (/> 1 ) x > x u 

(4) 

where r is the nominal temporal peak width of the 
pulse. The group velocity v g is given by 
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where the final limit applies when T 2 * ^> r. As 
the expression indicates, v g can be greater than c 
or even negative, depending on T 2 *, the width r 
of the pulse peak and g the atom-field coupling 



parameter. The expressions 
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gxo/2c and 



0i = r 2 g(xi — xo)/2c apply in the long T 2 * limit. 
The corresponding probability amplitudes for the 
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atoms in the medium are given by 
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C2(x, t) = — tanh— ( 
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Pulse and amplitude solutions were previously 
obtained by Andreev 19 for a fully coherent finite 
length medium in the absence of Doppler broaden- 
ing. The effects of Doppler broadening have been 
considered for infinitely extended media, 20 and an 
adiabatic EIT-based example has been given by 
Milonni. 21 

We define the pulse advance time as the amount 
by which the peak of a pulse moving through the 
gain medium precedes the peak of a pulse moving 
the same distance through vacuum. In our case we 
find 

L L L 2 _ 
Tadv — gr , (7) 

C Vg 2c 

where the group velocity v g in the medium is given 
in Eq. ([5]). Obviously, the pulse advance time can 
be substantial when the group velocity is negative. 
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Fig. 2. Snapshots of the analytic solutions given in 
Eq. (|4]). The horizontal axis is x in units of cr. The 
shaded region indicates the location of the medium 
which is taken to be 2cr in length. The parameters 
of the pulse and medium are described in the text. 
For an EIT comparison, see Fig. 15 in Ref. 21. 

We plot temporal snapshots of Eq. (j4]) in Fig. O 
where one can clearly see that the pulse entering 
the medium is advanced - in the third frame an 
outgoing peak leaves the medium even before the 
input peak enters it. In terms of initial temporal 
peak widths the peak advance is r a d v ~ 2.66r for 
a sample chosen to be Rb vapor, using the D2 line 
with these parameters: T 2 * = 0.733 ns, g = 266 



ns -2 and r = 0.1 ns, all of which appear to be 
within experimental reach. These parameters give 
an inverse Beer's length of a = y / 7r /2gT 2 * / ' c « 8.15 
cm -1 and a negative group velocity in the medium 
of Vg/c « —3.27. The medium was taken approx- 
imately 50 gain lengths deep, which corresponds 
to about 6 cm. Clearly Eq. Q predicts that the 
pulse advance time can be made arbitrarily large 
by simply increasing the length of the medium. 

4. Edge-Effect Instabilities 

After interacting with the medium, the analytic 
pulse solutions given in Eq. (|4|) return the atoms 
to their excited state, which is easily seen in Eq. 
((6]) where \c2(x,t — > ±00) | = 1. This is a fully 
coherent but unstable situation. We cannot realis- 
tically assume infinitely long pulse wings (as in the 
analytic sech functions). Since the studies of Som- 
merfeld and Brillouin, 22 the effects of finite cutoffs 
at the leading and also trailing edges of the pulse 
have been recognized as important to take into ac- 
count. Instabilities will always arise in gain me- 
dia from such effects. 23, 24 These instabilities can 
be understood via the two-level Area theorem. 14 If 
one defines the Area of the pulse to be 
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Q(x, t)dt, 



(8) 



then the well known Area theorem for an inhomo- 
geneously broadened gain medium (i.e. T 2 * <C r) 
is 

dO(x) a 



■ sm0(x), 



(9) 



dx 2 

where a = 7rgF(0)/c is the inverse Beers length. 
From Eq. ([9]) we can clearly see that pulses of Area 
2tt are unstable. The analytic pulse solutions have 
Area = 27r, so only the perfect infinitely long ana- 
lytic solutions are stable. Any slight variation, such 
as sharply cut off leading and trailing edges, will 
cause the pulse shape to change as the Area tends 
towards the stable values of tt or 3tt. 

We tested the sensitivity of the semiclassical 
equations to the cutoff instability by generating 
additional numerical solutions to the same set of 
equations, with the same medium parameters as 
in the previous section. We used a short enough 
pulse width to justify ignoring significant medium 
relaxation, but we gave the pulse weak but sharp 
leading and trailing cutoffs as shown in Fig. [3l In 
the previous section the analytic formulas predict 
a pulse advance time of about 2.5 peak widths for 
the parameters chosen, so we inserted cutoffs 10 
pulse widths from the pulse peak. 
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Fig. 3. The left frame shows the sech pulse used 
in numerical simulations. The right frame zooms 
in on the front edge of the pulse showing a sharp 
front edge 10 units from the pulse peak. The hor- 
izontal and vertical axes are the same as those in 
the analytic and numerical solution figures. 



Snapshots of the numerical solution are shown in 
Fig. [U and in frame 3 one sees an advanced peak 
leaving the medium before the input peak enters, 
along with a back-moving pulse (negative group 
velocity) in the medium, just as in the analytic so- 
lution. The output pulse is relatively cleanly sepa- 
rated from background, and is essentially identical 
to that for the analytic solution, in both its ex- 
tent of advance and in its Area which is 2tt. The 
negative pulse envelope, a kind of pulse "ringing" 
that trails the advanced pulse, is radiation arising 
from the coherence stimulated in the medium by 
the leading and trailing edges of the pulse. 
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Fig. 4. Snapshots during evolution of the numerical 
pulse solution to Eqns. ([2]) and ([3]) with the same 
parameters as for Fig. [2j but for a pulse initially 
entirely outside the medium with abrupt edges at 
±10 units from the peak. The pulse peak advances 
~ 2.5 units forward, which is still well behind the 
front edge of the pulse, which is marked by the 
arrows. 

The Area of the pulse calculated from the nu- 
merical solution is plotted in Fig. \5\ As the pulse 
enters the medium, its Area quickly changes from 



the unstable 2tt to the stable Area tt. However 
the Area does not settle to tt but rather oscillates 
slightly about that point. This is due to the fact 
that inhomogenous broadening does not dominate 
(cf. rephasing described by Burnham and Chiao 25 ). 
This Area change can also be seen in Fig. 0] where 
the analytic solution with Area 2tt is followed by 
the pulse ringing with Area oscillating about — tt. 

Thus, in spite of a cutoff instability, we have 
shown that the analytic pulse formulas in Eq. (j4j) 
and the analytic pulse advance time in Eq. ([7j) are 
still valid as long as the peak advance does not 
get ahead of the sharp front edge. After the pulse 
leaves the medium the instabilities appear as pulse 
ringing trailing the advanced pulse. 




Fig. 5. Time integrated area of the pulse from the 
numerical solution shown in Fig. [U The horizontal 
axis is position in units of cr and the vertical axis 
is the area of the pulse, as given in Eq. (|8|) in units 
of TT. 



5. Effect of Spontaneous Instabilities 

In addition to the cutoff instabilities mentioned 
in the previous section, purely quantum processes 
provide a different source for instability. Spon- 
taneous emission is the most fundamental, but 
single-atom spontaneous lifetimes are rather long, 
say 10's to 100's of nanoseconds, which we have so 
far overlooked. However, a more important quan- 
tum process is superfluorescence (SF), i.e., collec- 
tive spontaneous emission, for which the excited 
state lifetime is dramatically reduced, perhaps into 
the regime where it can't be ignored. Specifically, 
a superfluorescent pulse will emerge after a short 
delay time given roughly by td ~ 7tto/(47VA 2 L) 
for a sample with Fresnel number F ~ 1. Here TV 
is the density of atoms in the sample, To is the 
single atom radiative decay time due to sponta- 
neous emission, A is the wavelength of the emitted 
radiation, and L is the length of the sample. 26 ' 27 
This expression decreases with increasing length 



L, as should be expected, whereas the pulse ad- 
vance time given in ([7j) increases. Since successful 
observation of peak advance can be realized only 
if the SF pulse develops after pulse advance has 
occurred, this presents a conflict that must be ac- 
commodated. 

To examine this conflict we have used the results 
of Polder et al 28 who calculated the average delay 
time (average time for an SF pulse to emerge) to 
be given by 



(td) 



TTTqV 

4N a \ 2 L 



(ln(2irN a )Y 



3c 

4gL 
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where L = (27riVA) x / 2 . We plot both expres- 
sions (0 and JTDJ) in Fig. [6] as functions of medium 
length, in units of peak width cr. The result is pos- 
itive, in the sense that there is a range of medium 
lengths up to about L ~ 5cr for which the SF time 
is longer than the advance time. For these times we 
can expect to see pulse advances up to about 6.5 
peak widths. To ensure that the SF pulse delay is 
substantially longer than the pulse advance time, 
shorter lengths may be preferred, as we have cho- 
sen in our examples. 
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Fig. 6. Plot of the delay time, (td) of the superflu- 
orescent pulse given in Eq. ([TO]) and the total pulse 
transit time given in Eq. ([7j). The horizontal axis 
is the length of the medium in units of cr and the 
vertical axis is time in units of r. The dashed line 
is (td) and the solid line is r a d v - Each data point is 
an average delay time from 20 numerical solutions, 
showing excellent agreement with the theoretical 
delay time in Eq. ([TO]). Pulse advance is predicted 
for a medium of up to 5cr, giving pulse advances 
of up to 6.5 r. 

Spontaneous processes can be included in our 
semiclassical treatment by adding an incoherent 
stochastic trigger to the MS equations. We have 



done this via zero-average fluctuating dipole mo- 
ments as a surrogate for the vacuum fluctuations 
that are the fundamental basis for spontaneous 
emission. This modification is dealt with in detail 
in the Appendix. We then tested for pulse advance 
in the presence of both types of instability by solv- 
ing numerically the same MS system, equations ([2]) 
and ([3]), in this instance with edges cut off, and 
also with random initial dipole moments (see Eqn. 
([TT]) ). We have retained the same finite length 2tt 
sech input pulse described in Sec. IV. We used a 
medium 2cr in length (6 cm for a 100 ps pulse) 
just as in the previous examples. Actually all of 
the cases illustrated in this paper employ param- 
eters chosen in advance with guidance from the 
Appendix so as to satisfy the limitation that SF 
places on the length of the medium. The solution 
is plotted in Fig. and we see almost no differ- 
ence between that and Fig. [H where we assumed 
a perfect initial inversion. The pulse advance in 
this example is about 2.5 pulse widths in agree- 
ment with Eq. (|T|h and the advanced pulse exits 
the medium before a subsequent "ringing" pulse. 
The almost identical features of the solutions in 
Figs. [4] and [3 indicate that SF plays practically no 
role under the conditions chosen. 
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Fig. 7. Snapshots during evolution of the numerical 
pulse solution to Eqns. ([2]) and ([3]) with the same 
parameters as for Fig. [2j but with initial proba- 
bility amplitudes as given in Eq. ([TT]) . The pulse 
is initially entirely outside the medium with sharp 
edges at ± 10 units from the peak, with the front 
edge indicated by the arrows. The following peak 
advance is 2.5 units forward, in agreement with 
Eq. 0. This advanced pulse is also ahead of the 
superfluorescent pulse shown in Fig. 
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6. Summary 

We have examined theoretically the propagation 
of sech pulses in an idealized gain medium treated 
as inverted two- level atoms. The pulses are taken 
short enough to ignore collisional and radiative life- 
times for the atoms. Our motivation has been to 
test the feasibility of achieving substantial peak 
advance (so-called "fast light") in transmission 
through the medium, and here substantial means 
one full peak width or more. Fully coherent prop- 
agation of a resonant pulse in weakly Doppler- 
broadened Rb vapor was used throughout as the 
basis for our modelling. This is an example for 
which resonant input pulses are readily available 
experimentally with durations much shorter than 
the 26 ns spontaneous lifetime associated with the 
D2 transition. 

We exhibited an analytic solution to the non- 
linearly coupled semiclassical Maxwell- Schrodinger 
equations in terms of a pulse with 2tt Area, and 
showed in Sec. Ill that this solution is a "fast light" 
example with a peak advance that can theoretically 
be made arbitrarily large by adjusting the length 
of the gain medium. However we then included 
pulse cutoffs to introduce edges (which can travel 
no faster than the vacuum velocity). Nevertheless, 
a numerical solution carried out under the condi- 
tions employed for the analytic sech pulse showed 
that when these deviations are introduced, the fast 
light result is still obtained, and with nearly the 
same peak advance. A train of irregular pulses was 
induced by the edges, but they followed and did 
not interfere with the advanced peak. Finally, with 
the aid of modelling described in the Appendix, we 
took into account purely quantum effects. We used 
stochastic dipole fluctuations to induce collective 
spontaneous emission throughout the medium, i.e., 
superfluorescence, and found that it in our chosen 
regime it does not alter the peak advance in a sig- 
nificant way. 

In summary, we have extended the results ob- 
tained in previous thereortical examinations of fast 
light in fully coherent gain media. The combina- 
tion of results presented here suggests that pre- 
vious conclusions were unnecessarily pessimistic. 
Since the effect of superfluorescence is proportional 
to the number of participating atoms, it is desir- 
able to restrict the length of the medium, and in 
the Rb vapor example used, we could show that 
even in the presence of SF a peak advance of over 
one peak width could be arranged for transmis- 
sion of a 100 ps sech pulse in a 6 cm (50 absorp- 
tion depths) medium. Thus we predict that a very 



short input pulse can provide access to fast light 
using pulse bandwidths several orders of magni- 
tude greater than has been utilized to date, or pro- 
jected to be possible via EIT. 
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A. Maxwell- Schrodinger Model of Super- 
fluorescence 

To confirm that our conclusions in Sec. V satisfac- 
torily incorporate quantum fluctuations, we have 
compared results from our modeling with estab- 
lished theoretical studies of superfluorescence (SF) 
that have been backed up by experimental tests. 
The MS system that we are using assumes a clas- 
sical field. However SF can be extended to the MS 
system by appropriately modeling the quantum 
fluctuations with zero-average fluctuating dipole 
moments in the inverted medium. 28 30 Instead of 
perfect inversion, we assume the medium is ini- 
tially prepared as 

Cl (x^ = -oo) = sin(^)e^) 

c 2 {x u t -oo) cos ( — - — 1 , (11) 

where 0q{xi) is a randomly chosen very small ini- 
tial tipping angle of the Bloch vector, and (j){xi) 
is a random phase, for each point X{. For simplic- 
ity we assign <\>{xi) = or ir randomly to give a 
zero average complex dipole moment, ^(cicj) =0, 
when averaged over the length of the medium. Ar- 
bitrary phases in our numerical solutions do not 
change the results significantly, in agreement with 
previous studies. 29 

The modeling procedure assigns to the tipping 
angle 6q{xi) a very small randomly chosen num- 
ber with gaussian statistics at each point in the 
medium xi. The average value of the tipping an- 
gle is given by (0q) = with variance (0%) — 
(6o) 2 = 1/N a , where N a is the number of atoms 
and the average is over the length of the medium. 
We assume an initial density of atoms of N = 
8 x 10 10 cm -3 confined to a cylindrical atomic va- 
por cell with Fresnel number F = a /XL = 1. Thus 
the number of atoms is N a = NaL = NXL 2 ~ 10 8 , 
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for a 6 cm vapor cell. Since SF is a purely quan- 
tum mechanical effect, theoretical studies have de- 
rived these statistics for the initial tipping an- 
gle to ensure that the semi-classical MS system 
gives the same superfluorescent delay time (de- 
fined below) as purely quantum mechanical deriva- 
tions. 28 30 Subsequent experimental measurements 
confirmed these theoretical results. 16, 17 

Polder et al. 28 define the superfluorescent delay 
time, td, to be the time for the tipping angle in 
Eq. (jTTJ) at the output face of the medium, x = x\, 
to grow to become equal to 1. The SF delay time 
has alternative definitions, and has been derived 
in a variety of other ways, 26 ' 29 and the statistics 
have been studied. 30 The delay time specifically 
refers to the time it takes for the SF pulse to be 
emitted. However it can also be regarded as a decay 
time since, as the atoms emit the SF pulse, they 
subsequently decay to the ground state. 

We use the expression ([TP]) of Polder et al. for 
average delay time. To confirm that our numerical 
solutions conform to this model for superfluores- 
cence, we have averaged the delay times observed 
in 20 "experiments", i.e., for a given length L we 
made 20 different numerical solutions of Eqns. ([2]) 
and ([3]) using as an initial condition Eq. (fTTj) . but 
without any input pulse. An example of one of the 
solutions with a delay time of about (td) is shown 
in Fig. [HI This procedure was repeated for 9 dif- 
ferent lengths L, and the 9 averages are plotted 
as the data points shown in Fig. [6l One sees ex- 
cellent agreement between the average delay time 
from our numerical SF model and the expression 
([T0|) for the average delay time. 
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